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A theory is proposed for kinetic effects in isotropic Heisenberg antiferromagnets at temperatures 
above the Neel point. A metod based on the analysis of a set of Feynman diagrams for the kinetic 
coefficients is developed for studying the critical dynamics. The scaling behavior of the generalized 
coefficient of spin diffusion and relaxation constant in the paramagnetic phase is studied in terms of 
the approximation of coupling modes, ft is shown that the kinetic coefficients in an antiferromag- 
netic system are singular in the fluctuation region. The corresponding critical indices for diffusion 
and relaxation processes are calculated. The scaling dimensionality of the kinetic coefficients agrees 
with the predictions of dynamic scaling theory and a renormalization group analysis. The proposed 
theory can be used to study the momentum and frequency dependence of the kinetic parameters, 
and to determine the form of the scaling functions. The role of nonlocal correlations and spin-liquid 
effects in magnetic systems is briefly discussed. 
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I. INTRODUCTION 



In recent years there had been considerable activity, both experimen|tal-|and theoretical, in the fi^4i of critical 
phenomena of antiferromagnets with Heisenberg and RKKY interactionsH -0. This increase of interest)3-l3 was stim- 
ulated by new experiments in high Tc superconductors and heavy fermion compounds. In particular, the critical 
spin fluctuations are supposed to be responsible for NoiijJ'ermi-Liquid low temperature behavior of the heat capacity 



and resistivity in CcCuq-xAux and Cei-xLaxRu2Si2U^ in the viciijite-of the quantum critical point. Besides, the 
Fermi-type Resonating Valence Bond excitations were introduced bjij,t2l in the framework of the Heiseahexg model 
to describe unusual magnetic properties of high Tc cupratesEj and Ce - based heavy fermion compoundsO,E3. In this 
case the critical spin fluctuations were shown to be important in the spin liquid formation mechanism. However, the 
scaling behavior of kinetic coefficients in the presence of spin - liquid fluctuations can be differ signiflcantly from the 
predictions of the dynamic scaling theoryllj. 

In this paper we present a microscopic description of the scaling behavior of spin diffusion and relaxation coefficients 
in isotropic Heisenberg antiferromagnet in the critical region fk'^f^ the Neel temperature. The scaling dimension of 
kinetic coefficients was predicted by Halperin and Hohenbergci-llj. They proposed the dynamic scaling hypothesis 
based on the assumption that the critical exponents of kinetic coefficients are the same on both sides of the phase 
transition point. The naicroacopic treatment of the spin diffusion in the paramagnetic phase of ferromagnet was 
proposed by S.V.MaleevE^I -EJ. Within the framework of this approach the approximations which are necessary for 
validity of the dynamic scaling hypothesis were established and momentum and frequency dependence of the spin 
diffusion coefficient was analyzed. We will follow Maleev's approach to derive the scaling behavior of kinetic coefficients 
for Heisenberg antiferromagnet. 

The qualitativ&^ijcl|Ure of the critical phenomena expressed in terms of wave number k and inverse coherence length 
is well knownliSEl. There are two regions: " hydrodynamic" regime deals with the longitudinal fluctuations <C 1 
of the order parameter - the staggered magnetization N ~ Ni — N2 and "critical" one considers the fluctuations with 

3> 1 (the momentum q = |k — Q| corresponds to small deviation of the momentum from the AF vector). However, 
there is an additional conserved quantity in AFM which is the total magnetization M = Mi + M2. The fluctuations 
of this parameter do not lead to any singularities in spin correlators. Nevertheless, there is a region of longitudinal 
fluctuations of Af fc^ <C 1 which is also can be called as " hydrodynamic" . The goal of our work is to investigate the 
spin correlation functions behavior in the paramagnetic phase of AF and to establish the relations between the kinetic 
coefficients in two fluctuation regions of the phase diagram. 

The dynamic of magnetization fluctuations is described by the van Hove microscopic diffusion equation in the long 
wave limit 

^ = i^oV^M, (1) 
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where Dq is the spin diffusion coefficient. This equation is derived from the conservation low for the total magnetic 
moment, since the operator M commutes with the Hamiltonian. The fluctuation behavior of nonconserved order 
parametet-is radically different. According to the general postulates of the thermodynamics of weakly nonequilibrium 
processesc3 the rate of change N is proportional to the conjugate thermodynamical force 

here Go - is the static susceptibility, and the kinetic coefficient Tq satisfies the condition Fq > 0. The relaxation 
function determined by the equation (]l|) is nearly uniform, therefore the gradient terms for relaxation processes are 
omitted. In contrast to relaxation the first nonzero contribution to the spin diffusion equation is proportional to q^. 
Although the average value of M is equal to zero on the both sides of the phase transition, the fluctuations of total 
magnetization around zero take place. However, unlike the FM case the diffusion mode is not the "critical" mode for 
AF transition. 

Thus, we are interested in the dynamical susceptibility of a cubic Heisenberg antifcrromagnet in zero magnetic field 
above the Neel temperature, 

<ij> 

and the dipolar interactionS is neglected. As is known the susceptibility is related to the retarded Green's function 
by the following equality: 

x{K^)^{9f^ofKis{K^), (4) 
where g is Lande splitting factor, is Bohr magneton and 

poo 

Xfs(k,c.) = z / dte-* < [Sm.Sl^m >, (5) 



1 y-g-.kR,j 



M =< So >, N Sq^^^ > . 

Using the equations ^ one can obtain an expression for the retarded Green functions in the fluctuation region in 
the following form: 

■i Dk'^ 

K^si^ ^ 0, = /C(k, u;) = Go(fc)-^-^ (6) 

i^^s(q=(k-Q)-.0,^) = £(q,^) = — — 1— ^ (7) 

Here Gq is the static spin correlation function. The equations (|^, ^ correspond to diffusion and relaxation regimes 
respectively. 

When T — \T — Tc\/Tc <C Gi (Gi is the Ginzburg number which characterizes the conditions when the Landau theory 
is valid) the fluctuations become important. The dynamic properties in this region can be described by Halperin - 
Hohenberg scaling hypothesis. According to this theory the equation for the dynamic susceptibility can be expressed 
in terms of the scaling function F: 

KisiKco) = Go{k)FiH,-^), (8) 

Thus, the dynamic critical exponent z, which characterizes the critical fluctuations energy scale u) ^ , is connected 
with the static critical exponent v k, 2/i deflned by the expression ^ ~ t~'^ . We assume the static behavior of 
susceptibility in the form Go(<;) ~ ^'^^^ . We also neglect the Fisher index 77 which determines the "anomalous 
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dimension" ta. This approximation is valid in 3D casdla. To derive the scahng properties of AF one have to introduce 
two types of the scaling functions: T\ and T2 

/C(k,c.)=Go(k)^l(fcC,-;^), /:(q,c^)-Go(q)^2(gC,;p^), (9) 

In principal, expressions |^) det ermin e the kinetic coefficients Dq and Fq in terms of coherence length. Moreover, 
the renormalization group analysealiljlj has shown that the kinetic coefficients are singular in the fluctuation region 

r— I The theory under consideration is a variant of the mode-mode coupling theory proposed by Kawasaki0 (see alsoE^- 
c3). We extend to AF systems the dynamic scaling method offered by Maleev for ferromagnets. As was mentioned 
above, the aim of our paper is to derive the scaling behavior and momentum and frequency dependence of T functions 
in the fluctuation region of AF. Assumptions which are crucial for the microscopic substantiations of dynamic scaling 
hypothesis will be also formulated. 



II. GENERALIZED KINETIC COEFFICIENTS 



We consider the dynamic susceptibility of cubic Heisenberg antiferromagnet in the fluctuation region above the 
Neel temperature. The equations |^) can be rewritten in more general form: 

KU^M = (10) 

There is a simple expression for the spin - diffusion coefficient in the hydrodynamical regime 

L>o = lim lim fc-2^(k, w)Go"^(k), (11) 

The generalized kinetic coefficient in the relaxation region is 7(k, w) = F(k, w). The different limits fc ^ 0,Ci; ^ in 
the expressions (^, ^ significantly depend og-the relations between k and w. These limits are completely analogous to 
those considered in the Fermi Liquid theorycil. In the treatment below we consider the the quasistatic limit of k ^ 0, 

|L^|/k2 ^ 0. 

As it was shown by Maleev, it is possible to determine the kinetic coefficients in terms of the Kubo functions for 
the operators S and 5*, where S = dS/dt, 

^(k,c.) = , (12) 

l + Goi(fc)$^5(k,c.)' 



and 



^AB (k, uj)^- [K^B (k, c^) - (k, 0)] , 



K^siK oj) = i dte'^' < [Ak(t), S-k(O)] > . 
Jo 

The formal expression (|lj) is exact, and it takes into account the nonlinear character of the relaxation forces. In the 
case of purely exchange interaction, S ^ k and 7 — <I>^^(k, w). The Kubo function in the denominator vanishes at 
k = 0, and we arrive at the expression 7 = ^55- However, in studying the dispersion of the kinetic coefficients, i.e. 
its dependence on fc, w, we cannot, generally speaking, neglect this function in the denominator. 

Using the properties of spin operators, it is not difficult to obtain the relations between different retarded spin 
Green's functions K§g{\i,uj), K^^{'k,uj) and K^^(k,uj) in the paramagnetic phase. These relations are result of 
dispersion relationgc3: 

co'KisiKLo) = [Kf^ik,Lo)~Kf^{k,0)] (13) 

In particular, it is seen from equation ( p^ ) that the function i^^^(k, w) possesses the same symmetry properties as 
spin - spin correlation function. 
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The spin current operators taken in the form 



''° Pl,P2 

can be connected with the Kubo functions for Matsubara frequencies w„ with the use of equation of motion: 

= --^ E[^(p +^)- np)Kp,s^+^si^ (15) 
p 

The spin - current correlators can be calculated diagrammatically and then the explicit equations ( |lO| - |l^ ) for the 
kinetic coefficients can be obtained. When deriving eq. (^4|), we have confined ourselves to the lowest terms of the 
expansion in powers of ka and have put Vl^(p) « pTca^a, where a ~ 1. Since, we are interested in the correction 
terms of order (fc^)^ and (fc^^! "^"^ consistently neglect the (fca)^ and (fca)^ corrections, because ^ ^ a. Thus the 
problem of derivation of the kinetic coefficients is reduced to the analysis of four - spin correlators with current 
vertices. The later problem can be solved by analytical continuation of the temperature Feynman's diagram to the 
upper half-plane of complex variable lo The Feynman diagrams corresponding to the spin-current correlation function 
are shown on figure 1. 

The "bare" poles of the Green's functions |^) lie on imaginary axis. Introducing fictitious quasi - particles 
"diffusons" and "relaxons" we obtain the self-consistent equations for kinetic coefficients and determine their scaling 
dimensions. 

The static susceptibility in the fluctuation region obeys the Ornstein - Zernike law 

Go(q) = K§siH,0) ^ ^j,, ^^^^l^^ , (16) 

(here A is a constant (A~l),r<Cl). There is no singularities of static susceptibility in the diffusion region fc^ <C 1, 
where Go ~ A/{2Tc). 

In Sec. Ill the relations between kinetic coefficients are established and scaling dimensions of kinetic coefficient 
are analyzed. Sec. IV is devoted to a momentum and frequency dependence of spin diffusion and spin relaxation 
constants. 



III. RELATIONS BETWEEN KINETIC COEFFICIENTS 



We begin with analyzing the set of diagrams which cannot be separated into two parts by cutting one line of 
interaction, and include these diagrams into the irreducible self-energy part of the spin - current correlation function. 
Using the definition of 7 and properties of the functions K we obtain an expression for generalized kinetic coefficient 
in terms of irreducible self - energy parts: 

1 p 7^f Jk,w)7(k,u;)7^f<J(k,^^) , 7^f Jk, tj)7(k, cj) 

7(k,c.) = -[Sf^(k,..) - Ef^(k, 0) + '^\ [ /^_;.,i.f . ][1 + Go-'-^-^^^TT^TTTT^TT^] (17) 



The equation ( p^ can be derived—both from analysis of diagrammatic series for the spin - current correlatorJl^, 
and directly from Larkin equationE£l,t3. Below we use the following notations for retarded spin Green's functions: 
K^^Ck, uj) = 7?.^(k, uj)Kgg{'k, oj), and - are the irreducible self - energies. To obtain the graphical representation 

for the irreducible part one have to replace_the "dressed" vertex in fig.l by the irreducible one. The estimation 



of TZ in the framework of mean - field theoryll^H] results in: 

7^- (fc^)(fca) < (fcO^ (18) 

Moreover, due to analytical properties of spin correlators, TZ^ ~ w. We assume that the expression for TZ contains 
the small parameter a/^ in a limit of small u not only in critical region and neglect this contribution below. Therefore, 
the generalized kinetic coefficient 7 can be rewritten in terms of irreducible self - energies only: 



7(k, CO) = ^i^U^, to) - I]f^(k, 0)) (19) 
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The general set of diagrams for the irreducible self - energy part in imaginary time can be classified by the 
number of intermediate states. First, we confine ourselves to the diagrams with two - particles intermediate states 
(fig-2) 

Sg(k, icu) = (^^^)V ^^(kA(^)(p, k, ZLU, ^e, ^{cu - e)))(kA(2)\p, k, le, - e),^UJ))x 

e p 

X KssiP,ie)KssOi-p,iuj-ie) (20) 

Replacing the sum over momentum p by the integration one should use the "cutoff' p ^ S,^^ as an upper limit of 
integration. This means that the integrations are restricted by the region close to the appropriate singularities (small 
p and small q for p ~ q + Q in the vicinity of AF vector Q). From the point of view of the energy dependence, the 
vertices A^^^ being the three - point functions are the analytical functions of all three frequencies uj and, in eachj-of 
the frequencies they have a cut along the real axis; they have no singularities in the physical region of variables uMB. 
Due to these properties it is possible to separate the vertex into a static part and dynamic contribution, disappearing 
in the limit uj Oc2l. We concentrate below on the static vertex. 

As it can be seen from the fig. 2(a), the static vertices correspond to the longitudinal processes of "diffuson" - 
"diffuson" or "relaxon"- relaxon" pair creation and annihilation, since in this case only the modes of the same origin 
couple together. On the other hand, the static Green function G(k) is not sensitive to the direction of the momentum. 
In the other words, the corresponding lines in the diagrams are not directed. For this reason, the processes of diffuson 
and relaxon scattering contain the same vertices as the processes |Qf|pair creation and annihilation. However, for such 
a scattering processes in the limit k the Ward identity holdalaj^il (fig. 3): 

A(2)(p,k,0)~aGoVap (21) 

As was mentioned, the region of integration in the second term is concentrated near the points p « Q. The contribution 
from the staggered magnetization fluctuations into the spin diffusion can be calculated by means of the substitution 
p = q + Q and, by use of the property dGg^/dp = dG^^ /d(\. 

Let us consider the diagram fig. 2(b). The external momentum of this diagram can be taken equal to Qaf- Therefore, 
we have two different coupling modes, "diffuson" and "relaxon". Thus, this diagram corresponds to the process of 
"diffuson" - "relaxon" pair creation. For this reason we can not use the Ward identity to analyze this contribution. 
However, the "bare" vertex (fig. 3) has the scaling dimension 

A(')(p,Q,o)^9y/ap^p 

Besides, in the magnetic Brillouin Zone the points and Q are equivalent. Taking into account the direction inde- 
pendence of interacting modes in fluctuating region one can assume that the processes of multiple scattering on the 
static fields do not influence the scaling dimension of the static vertex. |^ 

Using the properties of vertices it is not difficult to perform the summation over e and the analytical continuation 
in tiO. As a result we obtain the following expressions for kinetic coefficients 

/OO 7 O O 

^cothi^) ^(VGo i(p))2 [/m/C(p, e) — ImlC{p - k, e) + ImC{p, e) — ImC{p - k, e)] (22) 
OO p 

^cof/i(^)X^(VGo-i(p)Q)2 X [/m/C(p,£) — /m£(p-q,e) + /m£(p,e)— /m/C(p-q,e)] (23) 

There index stands for the processes with only two - particle intepaediate states taken into account. The equations 
generalize the corresponding equation ( P3|) are obtained by Maleevllj for ferromagnets. In that case it is possible to 
retain only the first term in expression (|22|), since one deals with a single mode regime. Equations ( |2^ , |2^ ) can be 
rewritten in a different form. Taking k = and q = and integrating by parts over e one can obtain the following 
result: 



^As was mentioned above, the index of anomalous dimension (Fisher critical exponent) is taken to be zero. 
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These expressions are in fact a generahzation of equations derived by Maleev0 for the case of two couphng modes. 
As it was mentioned above, the main contribution to the integrals (^ 23) appears from the vicinity of scaling 



functions singularities (||). Besides, the characteristic energies of critical fluctuations satisfy the condition uj* ^ Tc 
due to the "critical slowing down". Thus, we can replace the hyperbolic functions (p2| -p4|) by their arguments. 
Calculating integrals over the frequencies and extracting the scaling dimension of kinetic coefficients we obtain the 
relation between spin diffusion coefficient and the relaxation constant. 



= b,T,'a%^r'— + b^T^a'i^)— (26) 
a JJq a 1 

It is interesting to note that the expressions (^) can be derived by direct substitution of retarded Green functions 
(|,0) to the expressions (|22|j2|). After integration over the frequencies we obtain the equations on kinetic coefficients 



which contain only integration over momentum of the different powers of the static correlator Gq. As a result the 
first term is determined by the two - diffuson intermediate state and the second term corresponds to two - relaxon 
intermediate state. 



The integrals (M - 25) can be calculated in the same way. The equation which determines the relaxation constant 



in terms of spin diffusion coefficient can be written in following form: 



To = ci(-)— +C2(-) (27) 



a To a 

The coefficients 61,2, ci_2 ~ 1 in eqs. ( p6p^ ) depend on the shape of a scaling function and, strictly speaking, can 
not be calculated in the framework of method used. Solution of self - consistent system of algebraic equations leads 
to following scaling dimension of kinetic coefficient |^ 

Do/{na^)r^Tor^i^/a)'/' (28) 

Behavio r pf . kinetic coefficient agrees with predictions of dynamic scaling theoryEUll and with renormalization group 



analysial^E2l. Therefore, the kinetic coefficients of AF are singular in fluctuating region, and spin diffusion is determined 
by the intermediate relaxation processes. A correction in (|2^ ) corresponding to the influence of self - diffusion is 
proportional to SDo/Do ^ i£,/a)~'^ t^^^ ^ 1- Hence, these processes can be omitted in the critical region. 
Calculation of dynamic critical exponent z (see eq.||) results in z = 3/2. 

The simple physical assumptions explaining the spin diffusion and relaxation in the ffuctuation region are based 
on the following argument: in the vicinity of the phase transition point T ^ the regions of short range order 
are formed. The characteristic length scale of the ordered regions is ^. The excitations in these regions are the AF 
magnons with sound like dispersion law. Estimations of the spin diffusion coefficienLresults in Dq ^ ^'^/^diff, where 
tdiff ~ £,/c is a characteristic diffusion time and c ~ jg a "sound" velocitjiiJ, therefore Dq ^ S,^^^- Due to 

the dynamic scaling hypothesis which postulates the invariance of the dynamic scaling exponent z, the characteristic 
scaling behavior of relaxation constant is Fq ~ 

Although the kinetic coefficients are singular, the characteristic relaxation times of staggered magnetization and 
spin diffusion become infinite. This lusault guarantees the existence of macroscopic ordered regions, corresponding to 
the partial equilibrium of the systemE!]. 

It should be emphasized that we did not use any assumptions about the elementary excitations spectrum in ordered 
phase when deriving the expressions (p6[). Our study is based only on a concept of total magnetization conservation 
and order parameter non-conservation. 



Spin diffusion coefficient in ferromagnet is not a singular function. It behaves as Do/{Tca^) ~ (^/ffl) 



-1/2 
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IV. CORRECTIONS TO THE KINETIC COEFFICIENTS 



Here we study the kinetic coefficient as a functions of finite but small k and u. For this sake we use the relations 
between spin correlation functions jl^ ) and Kubo functions (0). It can be seen that the corrections due to momentum 
and frequency dependence of kinetic coefficient are determined both by the nonlinear character of the relaxation forces 
and the momentum and frequency dependence of irreducible self - energy parts. According to the estimation ( [T^ ) 
one can derive the frequency and momentum dependence of spin diffusion and relaxation constants only in terms of 
linear response theory. Therefore, we can neglect the nonlinear character of relaxation forces. 

First, we consider the static renormalizations of kinetic coefficients. It is not difficult to obtain the usual expansions 
for spin diffusion and relaxation constants in series of (fci^)^", (g^)^": 

Z?(2)(k,0) =Z?o(0,0)[l + a (fc^^^...] 



r(2)(q,o) = ro(o,o)[i + /3'(gO' + -] 

As is knownil, these expansion can be done due to the singularities of static spin correlation functions in points 
ki = — n^^^^, where n are integer. The coefficients a , f3 are determined by the shape of static correlation function. 

We analyze the energy dependence of kinetic coefficients. Using equation ( ]T9| ) we obtain the following expressions 
for the real and imaginary parts of the generalized coefficient 7(k, co): 

Rej{k,Lu) = ^ , Im-/{k,uj) = (29) 

Since Imj is the odd function of uj and Re^ is the even function of uj the first nonzero term in the energy expansion 
of kinetic coefficients is proportional to cj^. 

Let us introduce the effective generalized kinetic coefficient 7* by the definition: 

9 , 



^ (30) 



1 + Go i(k)^^i?eSf<j(k,c^)U=o 



The expression for the coefficient 7* is completely analogous to the definition of effective mass in the theory of 
quantum liquids. The renormalization constant Z is similar to the Fermi - liquid Z factor. 

1 



1 + Go ^k)— ^i?eI]f^(k,c^)|^=o 



d_ 

Calculating Z in fluctuation region results in the following expressions for renormalization constant: 



where the coefficients e' , (5 <C 1 also can be expressed in terms of the integrals of different powers of static correlator 
Go. Using the definition ( |30| ) for small, but nonzero frequencies one can obtain the expansions of real spin diffusion 
coefficient D* and relaxation constant F* |^: 

D('>{k,u;)=DoiO,0)[l + aikO' + al^iu;/Lj*)' + -] 

F(2)* (q, u;) = Fo (0, 0) [/3 + /?' (g^^ + p'^^ {oj/u;*)' + ■■■] (32) 

It should be noted that we do not intend to describe the kinetic coefficients properties in the region w ~ w*, 
~ The diffusion and relaxation picture in this region becomes very complicated and scarcely amenable to 



Tct"^ is characteristic fluctuations energy, 2 = 3/2 
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detailed analysis at the present time. For this reason we shall neglect the irregular paixections to kinetic coefficients 
resulting from the generation of poles chain and drift of the poles toward the new cuul3,tll. In the region of k, uj under 
consideration, all these singularities are small corrections, and therefore, the phenomena associated with them are no 
special interest. 

Now we turn to discussion of the role of many - particles (to > 2) intermediate states. As it was mentioned above, 
we are interested only in regular contributions of such diagrams to the kinetic coefficients: 



{kaf J2 -Y. ^^"^k, pi, p„OA<™^\k, pi, p,„)^(pi + ...p„ - k)x 



pi 



dei.. .d£^ImK§s (pi , £i ) . ..ImK§g (p,n , e,-, 



-5{£l + ...+£m-Uj), (33) 



Here the functions K correspond both to " diffusons" and " relaxons" . The integration over frequency is still restricted 
by the regions close to singularities of scaling function. In the case of m = 2 the expression ( ^ ) simply transforms 
into (||,^. 

The vertices A^™^ are symmetric fuactions of the momenta of intermediate particles. For k = the generalized 
Ward identities similar to thous of refJl3 hold for A*^™); therefore, A^™' can be expressed as a sum of jderivatives of the 
ordinary n - particle vertices of static scaling theory. A "dimensional" estimate exists for the latteiO: r„j ^ p3-m/2 
in the limit k — > 0. As a result, the diagrams with m-particle intermediate states in the channel of "diffusons" and 
"relaxons" creation do not change the scaling dimension of irreducible self - energies. As to the vertices behavior 
around the AF vector Q, there is a direct mapping of the properties of the diagrams with two - particle intermediate 
states onto those with many - particle intermediate states. Therefore, we can confine ourselves only by the diagrams 
with two - particle intermediate states. Another diagrams only renormalize the numerical constants, and the latter 
can not be determined in a framework, of our method. The energy dependence of vertices also can not change the 
scaling behavior of kinetic coefficiental3. 



CONCLUSION 



We investigated the scaling behavior of generalized kinetic coefficients in 3D isotropic Heisenberg antiferromagnets. 
In a framework of analysis based on the modified mode - mode coupling theory the critical exponents are derived. 
Necessary approximations in microscopic approach to satisfy the dynamic scaling hypothesis are formulated. In 
particular, it is shown that the main sequence of Feynman diagrams contain the processes with two - particles inter- 
mediate states. Besides, the vertices of static scaling theory are proved to give the main contribution to quasiparticles 
scattering processes. 

The momentum and frequency dependence of spin diffusion and relaxation coefficients is found in the "pole" 
approximation. A definition of the generalized effective coefficient is introduced in a full analogy with that used in 
the theory of quantum liquids for a definition of effective mass. Taking into account the renormalization processes 
responsible for multiple scattering of "diffusons" and "relaxons" we derived the asymptotics of scaling function in the 
region of small frequencies and momenta, oj <^ uj* , k,q <^ ■ 

Our investigation of the critical dynamic is based on the static and dynamic scaling hypotheses. We took into account 
only two coupling modes. One of them is considered to be responsible for fluctuations of the total magnetization and 
another deals with the fluctuations of order parameter - staggered magnetization. Existence of two modes is determined 
by the conservation's laws of the Heisenberg model. For this reason, our results depend only weakly on the specific 
properties of AF and we expect that all they are valid for any system with nonconserved order parameter when the 
additional conserved quantity exists. 

The scaling behavior of more complicated systems, e.g. Ce - based heavy fermion compounds with almost integer 
valence, can be significantly different from critical dynamic of Heisenberg magnets. In these systems an additiojaal 
coupling modes interacting with paramagnons can appear due to, e.g., the spin liquid correlations, proposed in Ref.Ej. 
An additional modes in Kondo lattices can be also connected with itinerant - magnetic excitations. Kinetic coefficients 
of such systems can be measured by the methods of of inelastic neutron scattering to verify a spin liquid correlations 
existence. n n n 

The methods developed in this paper can be useful for the systems close to the quantum critical pointl3,Eil -E3 when 
the temperatures of AF ordering are close to zero, or, even negative. Our method also should be valid for anisotropy 
magnets, antiferrimagnets and systems with dipolar interactions. 
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We think that the diagrammatic approach for inveatigation of kinetic coefficient in the vicinity of Neel point has 
some advantages in comparison with other methodalj-tJm,-| This approach can be apphed to the unrenormahzible 
Hamiltonians and to systems with nonlocal coupling modeala. 
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FIG. 1. Diagrammatic equation for spin - current correlator. 
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FIG. 2. Feynman diagrams for kinetic coefficients. Only two - particles intermediate states are taken into account. The 
diffusion mode is represented by the wavy lines. Dashed lines stand for the relaxation mode. Dots represent the vertex parts 
from static scaling theory. 






FIG. 3. Equations for two - particle dynamic vertex. 
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